Hadron-quark phase transition in asymmetric matter with dynamical quark masses 



G.Y. Shao/ M. Di Toro/^^,* ^ Liu,^'^ M. Colonna/ V. Greco/'^ Y.X. Liu,^-^ 
^ INFN-Laboratori Nazionali del Sud, Via S. Sofia 62, 1-95123 Catania, 



o 

(N 
D 



> 

O 



^Physics and Astronomy Dept., University of Catania, 
'IHEP, Chinese Academy of Sciences, Beijinq, China 



and S, 
Italy 



Plumari^' ^ 



Theoretical Physics Center for Scientific Facilities, 
Chinese Academy of Sciences, Beijing, China 
^Physics and Astronomy Dept., University of Catania 
^Department of Physics and State Key Laboratory of 
Nuclear Physics and Technology, Peking University, Beijing 100871, China 
''Center of Theoretical Nuclear Physics, 
National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 

The two-Equation of State (EoS) model is used to describe the hadron-quark phase transition in 
asymmetric matter formed at high density in heavy-ion collisions. For the quark phase, the three- 
flavor Nambu-Jona-Lasinio (NJL) effective theory is used to investigate the influence of dynamical 
quark mass effects on the phase transition. At variance to the MIT-Bag results, with fixed current 
quark masses, the main important effect of the chiral dynamics is the appearance of an End-Point 
for the coexistence zone. We show that a first order hadron-quark phase transition may take place in 
the region T C (50 — 80) MeV and pa C (2 — 4) po, which is possible to be probed in the new planned 
facilities, such as FAIR at GSTDarmstadt and NICA at JINR-Dubna. From isospin properties of the 
mixed phase some possible signals are suggested. The importance of chiral symmetry and dynamical 
quark mass on the hadron-quark phase transition is stressed. The difficulty of an exact location of 
Critical-End-Point comes from its appearance in a region of competition between chiral symmetry 
breaking and confinement, where our knowledge of effective QCD theories is still rather uncertain. 

PACS numbers: 12.38.Mh, 25.75.Nq 



I. INTRODUCTION 

The determination of the phase diagram of strongly in- 
teracting matter and the search for signals of the hadron- 
quark phase transition are challenges both in theory and 
experiment. Intensive studies on these fields have been 
developed in the last decades [1-15]. Most phase dia- 
grams have been derived from Monte Carlo calculations 
of Lattice QCD [1, 4, 8] or effective chiral models [5~ 
7, 15] with quark degree of freedom. In the chiral effec- 
tive model, one can describe well the line of chiral phase 
transition and the complicated phase diagram of color 
superconductivity [16-18], while the confinement phase 
transition for the gluon part can be investigated with 
the extended NJL model coupled to the Polyakov loop 
[9, 12-14, 19-22]) in the temperature and chemical po- 
tential plane. However a large uncertainty remains about 
the derivation of baryon density and energy density of the 
onset of the hadron-quark phase transition [14, 23]. 

Although there arc attempts to describe nuclear and 
quark matter in unified effective models [24-28], further 
investigations are needed to give more satisfying results, 
in particular at high baryon and isospin densities, of in- 
terest for the expected phase transition in heavy ion col- 
lisions and compact stars. 

An alternative approach to describe the hadron-quark 
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phase transition is based on a two-EoS model with the 
Gibbs criteria, which has been widely used to make pre- 
dictions on the phase transition in the interior of neutron 
stars (e.g., [29-35] ). 

We remark that, even in the two-EoS approach, only 
a few papers have studied the phase diagram of hadron- 
quark transitions at high baryon density in connection 
to the phenomenology of heavy-ion collision ten A GeV 
range (intermediate energies) [36-40]. In Ref. [36] the 
phase transition from hadron to quark matter has been 
firstly analyzed for isospin asymmetric matter. It should 
be noticed that recently increasing attention is paid to 
this aspect, and some observable effects are suggested 
to be seen in charged meson yield ratio and in the on- 
set of quark number scaling of the mcson/baryon ellip- 
tic flows [37, 38]. This provides us a new orientation to 
investigate the hadron-quark phase transition, and can 
stimulate some new relevant research in the field. Later, 
hyperons have been included in [39], but the calculated 
results show that strange baryons are not important due 
to the exiguous final population in the short time scale of 
a nucleus-nucleus collision. Furthermore, the Cooper pair 
effect of w, d quarks (two flavor superconductor, 2SC) has 
been considered in the quark phase [40], which reduces 
the symmetry energy difference between hadronic and 
quark phases. The most important conclusion emerging 
from these works is that the onset density of hadron- 
quark phase transition is smaller in asymmetric matter 
than that in symmetric matter, which is possibly reached 
through heavy-ion reactions at new planned facilities, 
such as FAIR at GSTDarmstadt and NICA at JINR- 
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Dubna, where heavy ion beams (even unstable, with large 
isospin asymmetry) will be available with good intensities 
in the 1-30 A GeV energy region. 

One drawback in all these calculations is that current 
mass (or massless) u, d quarks are taken for the quark 
phase, where the MIT-Bag like models are used. The ob- 
tained results are possibly reliable at high density, after 
the restoration of chiral symmetry due to the asymptotic 
freedom of QCD, but the chiral symmetry breaking at 
finite temperature and low densities is not accounted for. 
From p A pq to smaller densities, the dynamical masses 
of M, d quarks become larger and larger, and they al- 
most reach 200 MeV at p ~ po for some parameter sets. 
This means that non-perturbative effects become more 
and more important with the decreasing baryon density. 
Therefore, calculations of properties of the phase diagram 
at low density and finite temperature arc not fully con- 
sistent. As a matter of fact when the MIT-Bag model is 
taken for the quark part, we see that the T-p, T-/i phase 
diagrams highly depend on the values of the bag constant 
_B, which cannot be determined accurately [38]. 

In order to obtain more reliable theoretical results and 
predict possible observablcs in the planned experiments, 
we take the Nambu-Jona Lasinio (NJL) model to describe 
the quark phase with the interaction between quarks, 
where the chiral symmetry breaking and restoration are 
well described. 

It was also proposed in Ref. [41], that in high en- 
ergy heavy-ion collisions strange and antistrangc quarks 
can be produced by thermal excitation (with net 
strangeness being zero required by the conservation law 
of strangeness in strong interaction) and strangeness 
would be much more abundant in the quark component. 
Therefore we will also consider the thermal excitation of 



strange and antistrangc quark, but we need to keep the 
chemical potential of strange quark to be zero, ps = 0, 
before the beginning of hadronization in the expanding 
process to make sure of the net strangeness being zero. 
One mechanism of hadron production is quark recombi- 
nation [42-47] , which is out of the range of the discussion 
in this paper. 

The paper is organized as follows. In Section II, we de- 
scribe briefly the used effective Lagrangians and give the 
relevant formulas for the Rclativistic Mean Field (RMF) 
theory adopted for the hadron sector. In Section III, we 
present the calculated phase diagrams and compare them 
with those obtained in the MIT-Bag model. Besides, we 
present some discussions and conclusions. Finally, a sum- 
mary is given in Section IV. 



II. NUCLEAR MATTER, QUARK MATTER 
AND THE MIXED PHASE 

For the hadron phase, the non-linear Rclativistic Mean 
Field (RMF) approach is used, which provides an ex- 
cellent description of nuclear matter and finite nuclei. 
One can calibrate the hadronic equation of state at zero 
temperature and normal nuclear densities, and then ex- 
trapolate into the regime of finite density and temper- 
ature. Our parametrizations are also tuned to repro- 
duce collective flows and particle production at higher 
energies, where some hot and dense matter is probed, 
see [48] and refs. therein. The exchanged mesons include 
the isoscalar-scalar meson((T), isoscalar- vector meson(a;), 
isovector- vector meson(/ci) and isovector-scalar meson((5). 
The effective Lagrangian can be written as 
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where the antisymmetric tensors of vector mesons are 
given by 



In the RMF approach all effective meson fields can 
be expressed via their mean values, simply related to 
baryon and scalar nucleon densities. In this way only 
nucleon degress of freedom are left to describe dynamics 
and thermodynamics of the system [49, 51]. The nucleon 
chemical potential and effective mass in nuclear medium 
are 



(2) 



M* = M - g^a - ggT^^S, 



(3) 



where M is the free nucleon mass, = 1 for pro- 
ton and Tgjj = —1 for neutron, and p* is the effec- 
tive chemical potential which reduces to Fermi energy 

E 



p + M*' at zero temperature. The baryon 



and isospin chemical potentials in the hadron phase are 
defined as 



H _ fJ-p ~T~ Pr. 
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The energy density and pressure of nuclear matter at 
finite temperature can be derived as 
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where /i(fc) and /i(fc) are the fermion and antifermion 
distribution funetions for protons and neutrons (i = 
p, n): 



Mk) = 



l + exp{(i?*(fc)-K)m' 



l + exp{{E*{k) + t,*)/Ty 



(7) 



(8) 



The effective chemical potentials fi* are determined by 
the nuclcon densities 
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With the baryon density p = Pb — Pp + Pn and isospin 
density p^ ^ pp ~ pn, the asymmetric parameter can be 
defined as 
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As model Lagrangians for the hadron phase, NLp (the 
isovector scalar meson 6 being not included) and 
NLpS (with the S meson) will be used. The effective me- 
son couplings heve been chosen to reproduce good nuclear 
matter properties and even to represent a reasonable av- 
erage of the density dependence predicted by Relativis- 
tic Dirac-Brueckner-Hartree-Fock calculation [52, 53], 
See details in Appendix Al of Ref. [37, 38] and also 
Refs. [48, 49, 51]. A note about strangeness in the hadron 
phase can be found in Ref. [54]. 



For quark matter, we use the NJL model [57] to de- 
scribe the interaction between quarks which is respon- 
sible for the quark dynamics at intermediate energies. 
The NJL model describes well the mesons spectra and 
successfully explains the dynamics of spontaneous break- 
ing of chiral symmetry and its restoration at high densi- 
ties/ chemical potential [58-63]. The Lagrangian density 
in the three-flavor NJL model is taken as 



k=0 

-K[detf{q{l + 75)g) + det/(g(l - 75)9)]. (11) 



where q denotes the quark fields with three flavors, u, d, 
and s, and three colors; mg = diag(m„, m^, mg) in 
flavor space; Afe are the Gell-Mann matrices and G and 
K the four-point and six-point interacting constants, re- 
spectively. The four-point interaction term in the La- 
grangian keeps the SUv{3) x SUa{3) x Uv{l) x Ua{1) 
symmetry, while the 't Hooft six-point interaction term 
breaks the ^7^(1) symmetry. 

As an effective model, the NJL model is not rcnormal- 
izeable, so a cut-off A is implemented in 3-momentum 
space for divergent integration. We take the model pa- 
rameters: A = 603.2 MeV, GA^ = 1.835, KA'^ = 12.36, 
iTT-u,!! = 5-5 and = 140.7 MeV, determined by Re- 
hberg, Klevansky, and Hiifner in Ref. [64] by fitting 
/tt, M^, niK and to their experimental values. 

The thermodynamical potential in 3-flavor quark sys- 
tem in mean field approximation is 
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where C is a constant to be fixed by physics condi- Ei{k) = y'k'^ + M*^ is energy-momentum dispersion re- 
tions, A^c = 3 is the number of color degrees of freedom, lation of the quark flavor i, and pi is the corresponding 
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chemical potential. The dynamical quark masses and 
quark condensates are coupled by the following equa- 
tions: 



and 



M* =m,-AG^, + 2K<p,cl,k [i^j^k), (13) 



-2Nr 



d?k M* 



(l-n,(fc)-n,(fc)), (14) 



here ni[k) and ni{k) are Fermi-Dirac distribution func- 
tion of quark and antiquark 



ni{k) = 



1 



l + exp{(i?,;(fc)-M.)/n' 
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The pressure and energy density of the quark system 
are: 



i—u,d 
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If we define 




+ 2G {^.J + + 0,2) - 4K(j)u 4>d 4>s , 

(19) 

we will find -Be// acting exactly as an effective bag con- 
stant, just like the bag constant in the MIT-Bag model. 
The difference is that now -Be// depends on the interac- 
tion between different quarks as well as on density and 
temperature. The same holds true for the Fermi motion 
contribution. For all the integrations the cut-off is im- 
plemented, together with the irrelevant thcrmodynami- 
cal constant C of Eq.l2, as usual by requiring the energy 
density and pressure equal to zero in the vacuum, ie., 
eQ = pQ = at zero temperature and density. 

The number density of quark flavor i can be derived 
with the relation pi = —dflq/dp,i 

^{n,{k) ~ n,{k)) . (20) 

As we have mentioned in the Introduction, although the 
strangeness could be produced at finite temperature in 
heavy-ion collision, the number of net strange quark is 
zero before the hadronization takes place in the expand- 
ing process. So we must keep fis = in the calculation 



according to Eq. (20), but the strange quarks in loop 
diagrams do contribute to the u, d self energies via a six- 
point interaction, and then they can affect the dynamical 
mass of u, d quark according to Eq. (13). The corre- 
sponding definitions of density and chemical potential in 
quark phase are as follows 

pi^\{Pu + Pd). p^^pu-pd, (21) 



P-B = 2^^^u + Pd), P3^^{Pu-Pd)- (22) 

The asymmetry parameter for quark phase is defined by 
aQ = -4=3^^i^. (23) 

p2 Pu + Pd 

The above introduction is a separate description of the 
hadronic and the quark phase. The goal of the present 
work is to extend these studies to the hadron-quark phase 
transition, in particular to stress the effect of dynamical 
quark masses. 

Since there are two conserved quantity, baryon num- 
ber and isospin, during the phase transition, the Gibbs 
criteria (thermal, chemical and mechanical equilibrium) 
should be implemented for the mixed phase. General dis- 
cussion of phase transitions in multicomponent systems 
can be found in Rcf.[29]. 
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If we define x the fraction of quark matter in the mixed 
phase, the Gibbs conditions for the mixed phase are 



P"{PB,P3,T) = pQipB,Ps,T), 



(24) 



where pg = (1 ^ x)Pi3 + Xp'b the total baryon den- 
sity in the mixed phase, and the total isospin density is 
P3 = i^ — xjp^+xp's- With the initial condition of asym- 
metric parameter in heavy-ion collision, the global 
asymmetry parameter a for the mixed phase should be 



a 



Pb 

a"{x = 0) 



Xfl 
XP% 



(25) 



according to the charge conservation. 



III. NUMERICAL RESULTS AND DISCUSSION 

Before presenting the phase diagram of hadron-quark 
phase transitions with the NJL quark model, we firstly 
display some results with the MIT-Bag model in which 
free fermions are considered with the current mass to„ = 
TOd = 5.5 MeV. As for the bag constant, different values 
are used in literature, here we take (-Bmit)^^^ = 160 MeV 
and 190 MeV, respectively, for comparison. All that will 
be useful for general comments on the nature of the phase 
transition as well as for a better understanding of the 
dynamical mass effects in the NJL approach. 

We will pay more attention to the phase diagram of 
asymmetric matter. Indeed the onset density of the 
hadron-quark phase transition for asymmetric matter is 
smaller than in symmetric matter [36-40], then with the 
possibility to probe it in heavy-ion collision at the new 
planned facilities, for example FAIR at GSI-Darmstadt 
and NICA at JINR-Dubna, We focus the attention on ex- 
periments with neutron-rich stable heavy beams, where 
large intensities can be reached. Here the largest accessi- 
ble isospin asymmetry parameter a is just above 0.2 (e.g. 



we have a = QflTJ 



238 T 



U-|- U collision, see also the 



Table II in Ref.[39]). Therefore we will consider a = 0.2 
in our calculation, just as taken in Refs. [37-40]. Of 
course the use of more asymmetric unstable beams will 
enhance the isospin effects described here. 



A. Some results with the MIT-Bag quark EoS 

In Figs. 1 and 2 we plot the T — pB and T — pg phase 
diagrams . In each figure the solid curves are the phase- 
transition lines from nuclear matter to quark matter for 
(Smit)^/"* = 160 MeV (190 MeV), and the correspond- 
ing dash-dot curves are the transition lines to pure quark 



matter. NLp hadron EoS has been used in the calcula- 
tion. 

From the two figures, it easy to see that the decrease 
of bag constant reduces the onset density of quark mat- 
ter, and in general that the phase transition curve highly 
depends on the value of bag constant, which is one of 
the motivations of this paper to use the NJL model to 
investigate the phase transition. In the Fig. 2 we see a 
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FIG. 1: (Color on line) The T — ps plane of asymmetric mat- 
ter with the isospin ratio ct — 0.2 for (-Bmit)^''* ~ 160 MeV 
and 190 MeV. The region between the solid and dash-dot 
curve gives the binodal surface of the mixed phase, where 
the hadron and quark matter coexist. NLp parametrization 
is used for the hadron phase. 




800 1200 

FIG. 2: (Color on line) The T — ps plane of asymmetric 
matter with the isospin ratio a = 0.2 for (Bmit)^''^ = 160 MeV 
and 190 MeV. The curves are hke in the previous figure. NLp 
parametrization is used for the hadron phase. 



variation of the baryon chemical potential along the tran- 
sition from pure hadron to pure quark phase, that could 
be interpreted as an evidence of a continuous second or- 
der phase transition for asymmetric matter, see [36]. We 
like to comment this point that we will see also clearly 
in the following results with the NJL quark EoS. 
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The transition is of first order, with the presence of a 
coexistence mixed phase, in both cases of symmetric and 
asymmetric matter. For symmetric matter, with only 
one conserved charge pB , we have the expected disconti- 
nuities of thermodynamical quantities hke entropy, pres- 
sure, at fixed chemical potential ^b- The presence of 
two conserved charges, pB and p^, in asymmetric matter 
keeps the first order nature of the transition but changes 
some properties inside the mixed phase. Due to the 
presence of a new degree of freedom the interaction can 
choose the most energetically favored charge densities in 
each phase, at each relative concentration, in order to 
minimize the free energy of the system [29] . We can have 
an increase of pressure and chemical potential inside the 
mixed phase but the transition is still of first order since 
for each x fraction we have a discontinuity in the pB and 
Pa densities of the two phases in equilibrium. We note 
that this effect, directly related to the internal interaction 
in the two phases, leads to the important isospin distilla- 
tion mechanism (a*3 > q,^) inside the mixed phase, as 
shown later, of interest for possible experimental signals. 
Similar properties are present in the liquid-gas transition 
for dilute asymmetric nuclear matter, which is behind 
the multifragmentation processes at the Fermi energies 
[50, 51]. 

Before closing the discussion of MIT-Bag results we 
would like to add a few comments about the interpre- 
tation of the T — pb and T — pB phase diagrams of 
Figs. 1, 2. 

For the lower bag costant value {B'^/^ ~ 160 MeV) 
we clearly see a squeezing of the binodal surface (mixed 
phase region) in the T — pB plane at temperature T ~ 
60 MeV, corresponding to a baryon chemical potential 
Pb — 900 MeV, and a re-opening at higher temperature 
and smaller density (or chemical potential). The same 
is happening for the large bag constant, although less 
evident. 

We note that this effect is appearing just around a 
transition chemical potential equal to the nucleon effec- 
tive mass Pb — M* p. We have a simple interpretation: 

• The "opening" at low baryon density and high tem- 
perature. 

Here the transition pB is below the nucleon effec- 
tive mass (not much reduced at low densities) and 
the effective chemical potential p* even smaller. 
Thus in the hadron part Ep ^ > p* and the fermion 
distribution function of the nucleons will be rather 
small. Therefore pg, Eq. (9), will show a slow in- 
crease with p B , while the pressure, mostly of ther- 
mal nature, is increasing (even for the positive an- 
tifermion contribution). At variance in the quark 
sector Pb always much larger than the used cur- 
rent quark masses (5.5 MeV) and so we see a fast 
increase of to balance the hadron pressure. 

• The "re-opening" at high baryon density and low 
temperature. 



With the increasing of baryon density/chemical po- 
tential, the nucleon effective mass decreases and the 
transition pi < M* to pi > M* takes place. Thus 
the fermion distribution function of nucleons and 
Pg will show a fast increase with pB- The interac- 
tion part of the hadron pressure will also increase. 
In correspondence we will need a fast increase even 
of the p§ in order to keep the = balance. 

We remark that the first argument is not anymore 
valid in the NJL frame, since the quark masses at low 
density and finite temperature recover the much larger 
constituent quark values. Therefore we would expect im- 
portant qualitative differences in this low pB phase dia- 
gram region when we take into account the chiral mass 
dynamics, see Subsection C. 

B. Results with the NJL quark EoS 

In Figs. 3 and Fig. 4 we present the P ~ pb phase dia- 
gram, within the NLp — NJL two-EoS scheme, respec- 
tively for symmetric and asymmetric (a = 0.2) matter. 
For each temperature the mixed phase region is between 
the two sohd dots. 

Clearly the pressure is a constant in the mixed phase of 
symmetric matter. Fig. 3, just the same as with Maxwell 
construction. At variance in the asymmetric case. Fig. 4, 
we have a monotonous increase. It is interesting to note 
that pressure rising is faster in the first part of the mixed 
phase (more evident for the lower temperatures T=30, 60 
McV). This is due to the isospin distillation effect, i.e. a 
large pd ~ Pu asymmetry for reduced quark fractions, see 
next Subsection, which is increasing the quark pressure 
as we can expect from Eq. (17). 

From Figs. 3, 4, we can also see that the size of the 
mixed phase shrinks with temperature. Meanwhile the 
onset density becomes smaller, opening the possibility 
of probing the coexistence phase in heavy-ion collision 
at intermediate energies. This effect is further enhanced 
by the isospin asymmetry as we will discuss in the next 
Subsection. 



L Symmetry energy effects 

In order to study symmetry energy effects we have per- 
formed the hadron-quark transition calculations with the 
hadron EoSs NLp and NLpS which present rather differ- 
ent symmetry terms at high baryon density, stiffer when 
the 6 meson is included [49, 51]. 

In Fig. 5 we plot the (identical) T — p^ phase diagrams 
of symmetric matter. For asymmetric matter a = 0.2 in 
Fig. 6 we show the results for the parameter set NLp and 
in Fig. 7 for the parameter set NLpS. Now the mixed 
phase region, the binodal surface, is rather different. In 
general the onset density of the mixed phase in asymmet- 
ric matter is smaller than that in symmetric matter, simi- 



7 



700 



NLp-NJL a=0.0 




10 



FIG. 3: (Color on line) Pressure of symmetric matter as a 
function of baryon number density at different temperatures. 
The solid dots correspond to the limits of the mixed phase. 
NLp — NJL two-EoS calculation. 
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FIG. 4: (Color on line) Like previous figure, but for the 
asymmetric matter a = 0.2 case. 



lar to the results obtained by the MIT-Bag model [36-40] . 
In fact the mechanism is the same: for a given baryon 
density in isospin asymmetric matter we have a larger 
repulsion in the hadron phase since the symmetry term 
is less important in the quark sector for both MIT-Bag 
and NJL effective lagrangians. 

If the NLpS parameter set is used, the onset density 
will be further reduced, as shown in Fig. 7. This is nicely 
due to the fact that at high baryon density the symmetry 
energy of hadron matter with the parameter set NLpS 
is much larger than in the NLp case, see the following. 
All that also indicates that isospin effects can be very 
important for hadron-quark phase transition in heavy- 
ion collision in order to shed light on nuclear interactions 
in an hot and dense medium. 

In Figs. 5, 6 and 7 the curves corresponding to a given 
quark fraction x inside the mixed phase are also shown. 
We note that the isospin effects are mainly relevant in the 
initial part of the coexistence region. This is consistent 



with the interpretation of the pressure behavior in Fig. 4. 
We finally remark that such lower density zone can be 
reached even at relatively lower beam energies. 
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FIG. 5: (Color on line) Phase diagram in the T — pg plane 
for symmetric matter with the parameter set NLp. Of course 
the same curves are obtained with the parameter set NLpS. 
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FIG. 6: (Color on line) Phase diagram in the T ~ pg plane 
for asymmetric matter with the parameter set NLp. 



To understand the role of the symmetry term on the 
hadron-quark phase transition, we further investigate the 
symmetry energy Egy„i in the two phases, defined in gen- 
eral by [37, 40] 



(26) 



where = for hadron matter, a'^ for quark matter. 

From the hadron (RMF approach) and the quark 
(MIT-Bag, NJL) EoSs described before we can evaluate 
the symmetry energy as a function of the baryon den- 
sity at fixed temperature. We display the result in Fig. 8 
for both hadron and quark matter at T=0, 80 MeV. It is 
clear the symmetry energy of quark matter, only due to 
the Fermi kinetic term in both MIT and NJL schemes, is 
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much smaller than that of hadron matter, where also in- 
teraction isovector mesons contribute. Moreover in gen- 
eral we get a decrease of the symmetry energy at lower 
densities and higher temperatures due to a smaller con- 
tribution from the Fermi kinetic motion. 

The symmetry energy of hadron matter with the pa- 
rameter set NLpS is larger than that of the NLp case. 
As discussed before by comparing Fig. 6 and Fig. 7, this 
leads to a smaller onset density of the transition. It is in- 
teresting to note that this is a genuine relativistic effects 
since the scalar covariant nature of the isovector 6 meson 
contributes to increase the symmetry energy at high den- 
sities directly with a larger repulsion and indirectly via 
a splitting of the neutron/proton effective masses, with 
M* < M*, see details in refs. [49, 51]. 

We see that larger is the symmetry energy difference 
between hadron and quark phase, smaller is the onset 
density of the phase transition. Moreover a big variation 
of the symmetry term will strengthen the observable sig- 
nals of phase transition, which will be discussed later. 
Of course all these effects will be enhanced by the global 
isospin asymmetry of the system, suggesting the interest 
in experiments with very neutron-rich unstable beams. 

In Fig. 9 we plot the T — fig phase diagram for sym- 
metric and asymmetric matter with both the parameter 
sets NLp and NLpS. For the symmetric matter at a 
given temperature, /is is constant in the mixed phase, i.e. 
varying the quark matter fraction x, then we have only 
one transition line (the empty circles). At variance, like 
in the MIT-Bag calculation of Fig. 2, in the asymmetric 
case we have a monotonous ps increase with increasing 
X, so there are two curves in T — p^ plane presenting the 
start and end, respectively, of the transition from nuclear 
to quark matter. What is important is that the chemical 
potentials of the onset of the transition (x = 0) in asym- 
metric matter are always smaller than the corresponding 
ones of symmetric matter. This is more evident with the 
parameter set NLpS (dotted curve), which again shows 



the importance of symmetry energy effects. 
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FIG. 9: (Color on line) Phase diagram in the T — pg plane 
for symmetric and asymmetric matter with the parameter sets 
NLp and NLpS. Always NJL quark EoS. 



2. Isospin Distillation inside the mixed phase 

In Figs. 10 and 11, we show the variation of the isospin 
asymmetry parameters of hadron (a^) and quark (a'^) 
matter inside the mixed phase at various temperatures, 
with the global asymmetry a = 0.2. We clearly see the 
much larger values of a'^ when the quark phase starts 
forming, roughly for a x fraction between 0.0 and 0.4. Of 
course the effect disappears when the pure quark phase 
is reached (% = 1) where the global asymmetry is recov- 
ered just for the two charges (pb, pa) conservation. This 
is a nice Isospin Distillation effect ruled by the symme- 
try energy gap in the two phases, as confirmed by the 
enhancement in Fig. 11, where the more symmetry re- 
pulsive NLpS EoS is used for the hadron part. 
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We note that the isospin asymmetry of quark matter 
decreases with the enhancement of temperature. This 
is due to a general decrease of symmetry energy effects 
at higher temperatures (and lower densities), as we can 
clearly see in Fig. 9. 

The color pairing interaction at high density can re- 
duce the isospin distillation since it is energetically equiv- 
alent to the introduction of an effective symmetry repul- 
sion in the quark phase [38, 40]. However the isospin 
effects discussed before are still present. Moreover we 
can expect the pairing correlations to be less important 
with increasing temperature. Experiments focused to ob- 
serve isospin effects in the mixed phase, using neutron- 
rich heavy ion collisions at intermediate energies, appear 
very appealing. 



experiments. As suggested before in Ref. [37, 38], an in- 
version in the trend of the emission of neutron rich clus- 
ters, 7r~/7r+, A'°/A'+ yield ratios in high density regions, 
and an enhancement of the production of isospin-rich res- 
onances and subsequent decays may be found. Besides, 
there is a controversial point of view about the enhance- 
ment of the yield ratio A/p [65-69]. 

For instance, in [37] the reaction 238y_|_238y ^ isospin 
asymmetry a = 0.227) at 1 A GeV has been investigated 
in the RMF approach, and a rather exotic nuclear matter 
is formed with baryon density around 3 — Apo, tempera- 
ture 50 — 60 MeV, likely inside the estimated mixed phase 
region, especially with the parameter set NLpd for the 
hadron sector. 



T=0 

T=30 MeV 
T=60 MeV 
T=75 MeV 




FIG. 10: (Color on line) Asymmetry of quark matter and nu- 
clear matter inside the mixed phase at different temperatures. 
The parameter set NLp is used for the hadron phase. 



3. New effect of the dynamical quark masses: a 
Cntical-End-Pomt ? 

If we compare MIT-Bag (with fixed current quark 
masses) and NJL (with chiral restoration mechanism) re- 
sults for the mixed phase wc remark only one main differ- 
ence. As shown in Figs. 1. 2, phase diagrams inT — p^ 
and T — Pq planes are derived, using the MIT Bag model, 
with a critical temperature reached only at zero baryon 
density. At variance the binodal curves at high temper- 
ature and low density cannot be obtained with the NJL 
model and we see a narrowing of the coexistence region 
up to a kind of Critical — End — Point, Figs. 5, 6, 7 
and 9. This important result derives from two qualita- 
tive new features of the NJL effective theory: i) the quark 
masses variation due to the chiral restoration, ii) the de- 
pendence of the effective Bag-constant -Be//, Eq. (19), on 
temperature and baryon density. 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



'■. \ NLp8-NJL a=0.2 


T=0 




T=30 MeV 




T=60 MeV 
T=75 MeV 











0.2 



0.4 0.6 
X 



0.8 



1.0 



FIG. 11: (Color on line) Like in the previous Figure, but 
with the NLpS EoS for the hadron phase. 



This behavior of the quark isospin asymmetry inside 
the mixed phase of the hadron-phase transition will af- 
fect the following hadronization in the expansion stage, 
finally producing some observable signals in heavy-ion 
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FIG. 12: (Color on line) Dynamical masses of u, d quarks 
inside the mixed phase at different temperatures. The solid 
circles indicate the two density limits of the coexistence re- 
gion. Asymmetric matter with a = 0.227. The NLpS EoS is 
used for hadron phase. 



The two effects are jointly leading to a determination 
of an End-Point of the mixed phase: 
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FIG. 13: (Color on line) Effective proton masses in the 
hadron (upper) and quark (lower) matter inside the mixed 
phase at T = 60 MeV. The solid circles indicate the two 
density limits of the coexistence region. Asymmetric matter 
with Q = 0.227. The NLpS EoS is used for hadron phase. 
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FIG. 14: (Color on line) Like in the previous Figure, but at 
temperature T = 75 MeV. 
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FIG. 15: (Color on line) Baryon density dependence of the 
NJL Effective Bag "constant", {Beuf^^, at T = 60 MeV. 
The line ending with solid circles indicates coexistence region. 
Asymmetric matter with a = 0.227. The NLpS EoS is used 
for hadron phase. 



FIG. 16: (Color on line) Like in the previous Figure, but at 
temperature T = 75 MeV. 



• The chiral dynamics largely increases the quark 
masses at low densities and finite temperatures, see 
Fig. 12 where we plot the u, d effective masses in- 
side the mixed phase at different temperatures, for 
asymmetric matter. If we reach the limit of a nu- 
cleon effective mass in the hadron phase smaller 
than the corresponding combination of quark effec- 
tive masses, e.g. for protons M* < 2M* + M^, 
we would have unphysical solutions for the Gibbs 
conditions since from Eqs. (9), (20) we will get only 
Ps < results. This is confirmed by the Figs. 13, 
14, where we plot the proton effective mass in the 
hadron phase (solid line) and in the quark phase 
(dot-dashed line) inside the coexistence region at 
T = 60 MeV and T = 75 MeV. We see that at 
the higher temperature we are close to a crossing, 
indication of the lack of physical solution at higher 
temperatures, as already seen in the corresponding 
NLpS results of the previous Figs. 7 and 9. 

We have a final interesting comment about the 
quark mass splitting shown in Fig. 12 for asym- 
metric matter. The difference in the quark masses, 
with M^ < M* , is larger at the onset of the mixed 
phase, where the isospin distillation effects induces 
a larger difference in the two quark-antiquark con- 
densates. 

• Actually we clearly see that in the same (T, p) 
region we do not have solutions of the Gibbs condi- 
tions. This is due to the second qualitative new fea- 
ture of the NJL approach, the density and temper- 
ature dependence of the "Effective Bag Constant" , 
Eq. 19, also related to the dynamical quark mass 
variation. As a consequence at low densities and 
high temperatures, for small values of the Bag con- 
stant we cannot get mixed phase solutions since 
the hadron pressure (mostly thermal) cannot equi- 
librate the quark pressure in the coexistence zone. 
In Figs. 15, 16 we show the density dependence 
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of {BeffY^^ at T = 60 and T = 75 MeV (the 
solid circles give the mixed phase limits). For the 
T = 75 AleV case we see a sudden drop around the 
onset of the mixed phase, good indication of a lack 
of solution for larger temperatures. 

In conclusion the chiral dynamics seems to lead to a 
Critical — End — Point (CEP) of the first order hadron- 
quark transition, around T ~ 80 MeV and pB — 2po 
or /is ~ 900 MeV. Beyond this point, i.e. at higher 
temperatures and smaller baryon densities, we have to 
follow an approach with just one effective EoS able to 
describe both phases in order to check if we get a conti- 
nous transition. 

It is interesting to note that the CEP appears when 
the NJL quark masses have reached the constituent 
quark values. However we are aware that the NJL ap- 
proach, with only chiral dynamics, is not good at low 
densities and finite temperatures since confinement is not 
accounted for (in fact our B^ff goes to zero), just where 
we can expect that hadron degrees of freedom will start 
to play a role [63]. A further investigation including a 
confinement mechanism is certainly required in order to 
confirm the End-Point evaluation and to understand the 
related physics. Finally correlation effects, beyond the 
mean field approximation, would be also more important. 

In any case the fact that in the (T , pb) plane around 
the End Point we definitely get a derivative dT/dfi ^ 
(see Fig. 9) could be an indication that we are actually 
reaching a continous transition for chemical potentials 
and entropies, as suggested by the Clausius-Clapeyron 
Equation [38]. 

Finally two more remarks: 

• From Figs. 5, 6, 7 and 9, we see that the End- 
Point of the mixed phase is not depending on the 
global asymmetry of the system and on the sym- 
metry term of the used hadron interaction. This 
is consistent with the reduced effect of the symme- 
try energy at high temperatures and low densities. 
Fig. 8. 

• Isospin effects appear still relevant in the region 
just below the End-Point, Tc (50 - 80) MeV and 
Pb C (2 — 4) pq, accessible in the transient compres- 
sion stage of Heavy Ion Collisions at intermediate 
energies. 

The latter point is interesting even because in that 
phase diagram zone the NJL and MIT-Bag results are 
very similar, for similar bag constants Bmit — B^ff 
[38]. This makes more reliable the observed isospin ef- 
fects. 



IV. SUMMARY 

In this work, we investigate the hadron-quark phase 
transition in isospin asymmetric matter. We use a Two- 



EoS approach as in all previous calculations. The novelty 
of this study is to insert chiral restoration effects on the 
quark masses using the 3-flavor NJL model for interact- 
ing quark instead of the MIT Bag model. We obtain the 
binodal surface of a first order hadron-quark phase tran- 
sition in the region of p^ > 2po and temperature T less 
than about 80 MeV, available in heavy- ion experiments 
in the near future. 

The calculated results show that the onset density of 
the phase transition is lowered with the increasing asym- 
metry parameter a, property possible to be probed in 
the new planned facilities, for example, FAIR at GSI- 
Darmstadt and NICA at JINR-Dubna, with realistic 
asymmetries for stable beams, see [37]. 

The phase diagrams at high density are like the ones 
given by the MIT Bag, model with appropriate bag con- 
stant, since the chiral symmetry is restored and the quark 
masses approach the current values. But the dynamical 
quark mass becomes more and more important with the 
reduction of baryon number density and the increase of 
temperature, causing a reduction of the Effective Bag 
Constant. The resulting effect is more relevant just in 
the region T C (50 - 80) MeV and pg C {2- 4) po, avail- 
able with the new planed facilities, where some suggested 
important observables may be found as the signals of 
hadron-quark phase transition in the future. 

The most interesting effect of the use of a consistent 
chiral quark mass dynamics is the appearance of a kind of 
Critical-End Point for the first order transition, around 
T ~ 80 MeV and PB ^ 2po or pB ^ 900 MeV. Fur- 
thermore the obtained phase diagram exhibits a region 
with confinement but chiral restored symmetry, as ex- 
pected for the Quarkyonic matter [70]. 

At variance with the results with the MIT-Bag quark 
EoS, the phase diagram at lower density and higher tem- 
perature cannot be derived from the Gibbs conditions in 
a two-EOS model, which means that chiral symmetry is 
very important in the phase transition. This conclusion 
stimulates new efforts in the search of an unique EoS, 
with chiral dynamics and confinement mechanism, able 
to describe both phases in the region just above the sug- 
gested Critical-End-Point. 

In any case reliable results have been obtained in the 
interesting region Tc (50 - 80) MeV and C (2 - 4) po, 
that will undergo the test of experiments in a near future. 
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